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A Kitaev-Heisenberg- J2- Js model is proposed to describe the Mott-insulating layered iridates 
A2lr03 (A=Na,Li). The model is a combination of the Kitaev honeycomb model and the Heisenberg 
model with all three nearest neighbor couplings Ji, J2 and J3. A rich phase diagram is obtained 
at the classical level, including the experimentally suggested zigzag ordered phase; as well as the 
stripy phase, which extends from the Kitaev-Heisenberg limit to the J1-J2-J3 one. Combining the 
experimentally observed spin order with the optimal fitting to the uniform magnetic susceptibility 
data gives an estimate of possible parameter values, which in turn reaffirms the necessity of including 
both the Kitaev and farther neighbor couplings. 



Frustrated spin systems have long served as a relatively 
simple yet rich source of exotic phenomena such as spin 
liquids and unconventional order. The frustration may 
arise either geometrically on a lattice incompatible with 
the spin ordering, or dynamically from non-commuting 
competing terms in the Hamiltonian. The nearest neigh- 
bor S = 1/2 Heisenberg model on the kagome lattice is 
an instance of geometrical frustration that may even host 
a quantum spin liquid ground stat^. Bipartite lattices 
such as the honeycomb can still be geometrically frus- 
trated by including farther than nearest neighbor anti- 
ferromagnetic Heisenberg exchange, giving so called Ji- 
J2-J3 models. Such models on the honeycomb in par- 
ticular have seen a recent surge of worlPj^^, though a 
quantum spin liquid phase may require charge as well 
as spin fluctuations^H^. Breaking spin rotational sym- 
metry provides avenues for dynamical frustration, as in 
the Kitaev honeycomb model^^, a nearest neighbor Ising 
coupling of spin component set by a bond label 7 as in 
Fig.[l] This seemingly artificial model is exactly solvable 
with a spin liquid ground state exhibiting an emergent 
Majorana fermion with a Z2 gauge background. 

A recent and surprising addition to the experimentally 
relevant J1-J2-J3 models of frustrated spin systems, the 
Kitaev coupling has been recently proposed^ to oc- 
cur in the Mott insulatin^^ iridates A2lr03 (A=Na,Li), 
where the iridium ions are arranged in layers of 2D honey- 
comb lattices. Uniform suscep tibility and heat capacity 
studies on these materials^^'^ found Curie- Weiss temper- 
atures of —125 K for Na2lr03 and —33 K for Li2lr03, and 
a low magnetic ordering temperature of 15 K for both, 
suggesting strong frustration. A resonant x-ray scatter- 
ing measurement^^ on Na2lr03 found the ground state 
has antiferromagnetic order at wavevector M, suggested 
by a first principles calculation^^ to be a zigzag rather 
than a stripy configuration (see Fig.|2|. 

Strong spin-orbit coupling splits the iridium t2g states 
into a filled manifold and a half filled Kramer's doublet, 
an effective spin- 1/2 degree of freedom which need no 
longer respect the rotational symmetry. Thus the 90° 
angles of the Ir-O-Ir hopping path within the oxygen 
octahedra, together with d-orbital Hund's rule coupling 



and orbital interactions, are able to give the Kramer's 
doublet highly anisotropic exchanges of the Kitaev form. 
Higher order hopping paths, direct orbital overlaps, trigo- 
nal distortions and spin-orbit energy splittings within the 
iridium two electron propagator all contribute spin inter- 
actions other than the Kitaev term, primarily including 
antiferromagnetic Heisenberg exchange. 

Keeping only the nearest neighbor Heisenberg ex- 
change yields the Heisenberg-Kitaev modeJ^^lMisJ^ 
Eq. ([1]) with J2,J3 set to zero, which has been previ- 
ously used to describe the A2lr03 materialJ^^H^. The 
phase diagram^^ in the parameter < a < 1 con- 
sists of a Neel phase for the Heisenberg model at small 
< a < 0.4, the Kitaev spin liquid at large 0.8 < a < 1, 
and an intermediate antiferromagnetically ordered stripy 
phase (see Fig.|2|. The stripy configuration is the exact 
ground state at a = 0.5, solvable by means of a periodic 
site dependent spin rotation^^ which turns the Hamilto- 
nian into a Heisenberg ferromagnet in the rotated spins. 

Preserving J2 and J3 to produce the previously un- 
studied Kitaev-Heisenberg- J2-J3 model is important for 
two reasons. First, substantial J2 and J3 are likely to 
exist in the materials; density functional theory (DFT) 
calculations^^ for Na2lr03 found J2/J1 ~ 0.5, and a later 
tight binding fit of the DFT data including J3 found 
J2, J3 to be approximately equal^^. Second, the exper- 
imentally suggested zigzag ordered ground state^ can 
not be realized in a Kitaev-Heisenberg model alone. It 
is found that an antiferromagnetic J3 term is needed to 
stabilize the zigzag order. Moderate Kitaev and J2 cou- 
plings stabilize both zigzag and stripy orders. We will 
also show that in order to reproduce the experimentally 
measured uniform susceptibility x{T), the farther neigh- 
bor J2 and Js couplings as well as the Kitaev term are 
likely needed. 

The Kitaev-Heisenberg- J2-J3 Hamiltonian is 



H = J 
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where (ij), {{ij)) and {{{ij))) stand for the first, second 
and third nearest neighbor bonds, and jij is a nearest 
neighbor bond label, as illustrated in Fig.[l] The model 
interpolates between the J1-J2-J3 model at a = and 
the Kitaev model at a = 1, maintaining the second and 
third neighbor coupling strengths J2 and J3 in units of 
the nearest neighbor Heisenberg coupling strength. 

A recently proposed alternative model for Na2lr03 
based on ab initio calculations^^ takes the limit where 
trigonal distortion effects are stronger than spin-orbit 
coupling, finding a Hamiltonian with Ising anisotropy 
and no Kitaev term^^^. Putting this interesting scenario 
aside^^, we find that mild c = (1, 1, 1) uniaxial trigonal 
distortion is consistent with our approach. The effec- 
tive spin-1/2 Kramer's doublet remains well separated 
from the filled states. Its modified wavefunction creates 
anisotropics in the magnetic field coupling (^-factor ten- 
sor) and combines with the non-90° Ir-O-Ir hopping path 
to perturb Eq. (1), possibly enhancing both Kitaev and 
Heisenberg terms in addition to creating small Ising S^S^ 
and Ising-Kitaev S^^^ terms. Both modifications are 
expected from the observed anisotropy in single crystal 
Na2lr03 susceptibilit}'^^ and do not change our results. 
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FIG. 2: (a) Sample (J2, J3) slice of the classical phase dia- 
gram, with phases (I), (II) and (IV) represented in (b), (c) and 
(d) respectively. Region (III) contains various noncoUinear 
spiral configurations, (b) (I) Neel (c) (II) Zigzag, (d) (IV) 
Stripy. 
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FIG. 1: The honeycomb lattice, with Kitaev label 7 for near- 
est neighbor bonds, and including second and third neighbor 
bonds with Heisenberg couplings J2 and J3. 

Since there is ample evidenc^^^H^^ for magnetic order- 
ing in both Na2lr03 and Li2lr03 we will leave the cal- 
culation of the quantum phase diagram of Eq. for fu- 
ture work, instead turning to the magnetically ordered 
phases which may be studied by a purely classical anal- 
ysis. For each point (a, J3) in the three dimensional 
phase diagram we determined the magnetic ordering con- 
figurations using a quadratic (unconstrained) classical 
spin modeP^, which we diagonalized analytically in mo- 
mentum space. Since F = — F and M = —M these two 
wavevectors automatically give configurations of collinear 
unit-length normalized spins despite the absence of the 
unit-length constraint in the calculation, reaffirming the 
validity of the classical solution at these points. So- 
lutions at wavevector K or at generic incommensurate 
wavevectors correspond to noncollinear spiral configura- 
tions, which we label as a single phase. 

In order to discuss results on the classical phase dia- 
gram we introduce standard nomenclature from the lit- 
erature. For each ordering wavevector the phases are la- 
beled by a Roman numeral^ as follows. F: (I) Neel. M: 
(IV) strip'J^, and (II) zigzaj^ (or columnar^). All other 



wavevectors: (HI) spiral. Figure [s] displays six (72,^3) 
slices of the classical phase diagram at various fixed a. 

Quantum fluctuations modify the classical phase dia- 
gram in two ways. First, they create regions of quantum 
phases such as the plaquette valence bond solid or the 
Kitaev spin liquid; the former has been seen in the Ji- 
J2-J3 model^, while the latter appear j^^P^ at small J2, J3 
starting at a > 0.8. Second, they shift the boundaries 
between the magnetically ordered phases. (Qu antum fluc- 
tuations disfavor the spiral configurations^ ^ in favor of 
the collinear ordered phases, shrinking region (H I); th ey 
also favor the Neel state (I) over the other order^^^. 

The three dimensional phase diagram offers insights 
otherwise unavailable in its various limits. The stripy 
(IV) region in the J1-J2-J3 model at a = is in the 
same phase as the fluctuation-free exactly solvable point 
a = 0.5, J2 = J3 = which may be understood only 
within the Kitaev-Heisenberg model^^. As a increases, 
both the stripy and the zigzag phases grow substantially 
larger. The dynamic frustration by the Kitaev term and 
the geometric frustration by the J2 term have similar 
effects on the ordered phases, destabilizing Neel in favor 
of stripy and zigzag. 

It is worth reporting the direction of magnetic ordering 
in the various phases (excepting the special points a = 
and a = 1/2). The direction of the collinear magnetic 
ordering in both stripy and zigzag phases is constrained 
already at the classical level. For stripy order the 
spins lie along 5^, as was already determined by the 
spin rotation^^ solution of the J2 = -^3 = 0, a = 0.5 
Hamiltonian. For zigzag order we found that the spins 
are constrained to the S^S^ plane (see Fig.[2|. Thermal 
and quantum fluctuations ("order from disorder") force 
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the spins to lie along a cubic axis within the classically 
allowed space, in this case the and axes. Trigo- 
nal distortion gives other perturbations: for example for 
Mz stripy order it cants the spin axis from toward 
the distortion axis, and for the Neel phase the distortion 
axis may be an energy minimum or maximum within the 
Bloch sphere. A linear spin wave analysis found that di- 
rections closest to cubic axes are still preferred by quan- 
tum fluctuations. However, anisotropy in the real ma- 
terial likely overcomes all these effects to determine the 
ordering direction^^. 

Next we discuss the comparisons between experimen- 
tally measured susceptibilit}'^^ and exact diagonaliza- 
tion (ED), flrst describing each in turn. Uniform mag- 
netic susceptibility data for the sodium and lithium ma- 
terials at temperatures up to 300 K was taken from 
the most recent stud}^^"^, with the constant background 
removecF^. We used data from temperatures above 15 OK 
in order to avoid flnite size effects when comparing to 
ED. ED using the "fulldiag" ALPS moduli was per- 
formed keeping all eigenstates to enable comparison with 
high temperature data. The system diagonalized was an 
eight spin cluster, the unit cell of the a = 1/2 site de- 
pendent spin rotatioEpS, with periodic boundary condi- 
tions. As expected, the eight-spin ED, corresponding to 
a high temperature series expansion with eight-spin clus- 
ters, is reliable to far lower temperatures than the two- 
spin Curie- Weiss expression which only holds at T J. 
We found that ED flnite size effects for eight-spin clusters 
were only visible in the susceptibility at low temperatures 
T < J/2, well below J. The highest J values needed for 
good flts were below the 150 K data cutoff, self consis- 
tently affirming the reliability of the ED flts. 

For each parameter set (a, J2, J3) we diagonalized the 
system to generate a curve x(^)- The Hamiltonian 
Eq. ([1]) with a magnetic fleld coupling term has two pa- 
rameters in addition to (a, J2,-^3), namely the overall 
scale J and the magnetic fleld coupling g ja. Since the 
effective spin- 1/2 turns out to have the same ^-factor as 
an electron spin, we flx ^ = 2 and expect /i//iB to remain 
close to jii/jiiB = 1- For each (o^, J25^3) point the curve 
x{T) was flt to the experimental data by the two parame- 
ters J (corresponding to horizontal stretching) and fi/fiB 
(with {11/ (ib)'^ corresponding to vertical stretching). The 
resulting flt was evaluated by a "goodness function," the 
product of three Gaussian distributions, enforcing the 
following three conditions for a good flt. First, the mag- 
netic moment /i//iB found by the best flt must be close 
to 1, with a standard deviation of 0.15. This constraint 
on ji effectively constrained J as well. Second, the root- 
mean-square relative flt residual must be near zero with 
a standard deviation of 10~^. Third, the third neighbor 
coupling must be smaller or not much larger than the sec- 
ond neighbor coupling, J3 < J2, relaxed by a standard 
deviation of 0.2. The absolute (unsealed) value of this 
goodness function was used to produce the shading in 



Fig.[3j with darker shading corresponding to better flts. 

Given knowledge of the ground state magnetic order in 
Na2lr03 and Li2lr03, appropriate values for a, J2 and 
J3 are found by intersecting the darker shaded regions 
in Fig. [3] with the domain of the ordered phase. The 
estimated Na2lr03 and Li2lr03 parameters given either 
stripy or zigzag magnetic order are summarized in Table 
I. All material and order combinations yielded fltted val- 
ues of J in the range J ^ 60-150 K, with the likeliest 
values J ~ 100 K. The lithium material has less struc- 
tural distortion than the sodium material^^, suggesting a 
larger a, in agreement with the fltting results if they have 
the same magnetic order. For zigzag ordered Li2lr03 we 
flnd a ~ 0.7, i.e. Jk ^ 4-5 Ji with a numerical value 
of Jk ^ 130 K. Such a large Kitaev term relative to 
the other couplings suggests that the Kitaev spin liquid 
phase may be within experimental reach^^. In particu- 
lar, doping Li2lr03 may suppress its magnetic order to 
reveal characteristics of a doped Kitaev spin liquid^^. 

In conclusion, we propose the Kitaev-Heisenberg- J2- J3 
model, determining its ordered phases and further using 
ED flts of susceptibility measurements to demonstrate its 
applicability to Na2lr03 and Li2lr03. We flnd that the 
geometrical frustration due to J2, J3 and the dynamical 
frustration due to the Kitaev term both stabilize the 
same unconventional stripy and zigzag ordered ground 
states before the onset of the Kitaev spin liquid. We 
extract appropriate values for the spin couplings by 
flrst restricting to the experimentally observed magnetic 
order in the phase diagram, and then by requiring good 
fltting of the susceptibility x{T) by ED data. For zigzag 
ordered Li2lr03, a signiflcant Kitaev term Jk ^ 130 K, 
flve times larger than the nearest neighbor Heisenberg 
coupling, as well as substantial J2 and J3 couplings, are 
required for good agreement with experimental data. 
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ship Council and the Tsinghua Education Foundation in 
North America (Y.Z.Y.). 



TABLE I: Parameters for given M-wavevector order 





NasIrOs: 


Stripy (IV) 


a ^ 0.2-0.3, J2 < 0.5, J3 < 0.2 J ^ 110 K 


Zigzag (II) 


a ^ 0.4-0.6, J2, J3 > 0.4 J 100 K 




Li2lr03: 


Stripy (IV) 


a ^0.5, J2,J3<0.3 J;:^100K 


Zigzag (II) 


a ^0.7, J2,J3>0.4 J^90K 
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FIG. 3: (Color online.) Fixed a slices in ^J^^Jz) showing the magnetically ordered phases (I, II, III, IV) = (A^ee/, zigzag^ spiral, 
stripy) and shading corresponding to the ED x(T) fit goodness. Increasing the Kitaev term (i.e. increasing a) enlarges the 
extent of the zigzag and stripy phases, which occur at both small and large a. Fits to Na2lr03 are shaded in orange (with 
dotted contour lines) and fits to Li2lr03 are in blue (with dashed contour lines); darker shading corresponds to good fitting with 
m/ms ~ 1 and J2 > J3, while lighter shading corresponds to poor agreement. Given a magnetically ordered ground state for 
each of the materials, the range of allowed parameters is found by intersecting the darker shaded region with the magnetically 
ordered phase. 
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